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In recent work, the Schwarzschild interior of a black hole was investigated, incorporating quantum 
gravitational modifications due to loop quantum gravity. The central Schwarzschild singularity was 
shown to be replaced by a Nariai type universe. In this brief report we show that this interior 
solution is stable with respect to small perturbations, in contrast to the classical Nariai universe. 
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do: 

^^ , In a recent article [l[, the interior of a Schwarzschild black hole was analyzed in the framework of loop quantum 

^^ ' gravity (LQG - for reviews see [1, y, |3l)- In the analysis it was shown that loop quantum effects can cure the central 
^^ , singularity and replace it with a Nariai type universe with a Planck scale two-sphere radius. This would seem to imply 
bJO' that infalling matter would settle in a finite Planckian radius as opposed to collapsing to the classical singularity, 
^ I though a more rigorous analysis of the loop effects would be needed to answer such a question. However, concerns 
■^^ ' were raised in [5| that the Nariai solution of classical general relativity is unstable to certain perturbations with the 
instabilities leading to a Schwarzschild-de Sitter solution. 

In this brief report, we show that the Nariai solution in the interior is stable owing to the loop quantum effects 

in contrast to the classical situation. We perform a perturbative analysis about the Nariai solution using the loop 

^~^, quantum equations of motion. The linearized perturbation equations can be solved analytically and we show that 

pl^' the perturbations decay if evolved towards the center thereby establishing the stability of our solution. Note that 

I I the issue of stability we investigate is separate from the stability studied in [6|. There, the stability of the underlying 

discrete quantum theory is discussed in the context of various quantization ambiguities of the theory. The analysis of 

[6| is relevant to this work in that it supports the choice of S parameters that wc use in the effective Hamiltonian ^ . 

We briefly recall the equations describing the Schwarzschild interior in the LQG framework. For more details the 

f^ ' reader is referred to [3, Q. The interior has been studied in LQG in @, [l3, llll- For more details, the reader is referred 

•^ , to those articles. The interior metric is of the form 

g ; ds^^ ~N{tfde + g^.,{t) dx^ + gnn{t) d^^ (1) 

CO ' with N{t) being the freely specifiable lapse function and dQ" representing the unit two-sphere metric. The classical 
r — . Schwarzschild interior solution is given by 
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_ ^ • for t in the range t € [0, 2m] , where m is the mass of the black hole. Here, i = corresponds to the classical singularity 
►v> , and t = 2m being the horizon. Loop quantum gravity uses a separate set of variables than the metric ones. The 
; I ■ metric components are encoded in triad components Pb,Pc which are related as 



gxx = Y-^, gnn = \Pc\- (3) 

be I 

Note that the physical radius of the two-sphere is determined from the g^n or pc components. In particular, the 
classical singularity corresponds to a vanishing two-sphere radius where gnn = Pc = 0- The LQG variables also consist 
of connection components b, c which are conjugate to the triad components. We focus on the triad components for 
interpretation since they are directly related to the metric components. 

The Nariai solution [IJ, [l3] is a spherically symmetric solution of the classical field equations of general relativity 
where the two-sphere radius is assumed to be constant and the model is sourccd by a positive cosmological constant 
A with the possible inclusion of charge. For the case of no charge, the metric is given by 

ds^ = -dt^ + 4 cosh^(%/At)rfa;2 + ^dfl^ (4) 

A A ' 
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(5) 



the form of which is equivalent to ^ . The triad components are thus given by 

— cosh(vAt), Vr ~ -r ■ 
A ^ ^ A 

The classical Nariai solutions are quantum mechanically unstable. Such instabilities lead to perturbations in the two- 
sphere radius Pc and the solution decays into a Schwarzschild-de Sitter black hole solution (or Reissner-Nordstrom-de 
Sitter for the charged case) [ij, E, 13 • 

The loop quantum effects are derived from a modified effective Hamiltonian [l|, |T3| based on the phase-space 
variables 6, c,pb,Pc which is given by (note that we use the improved quantum Hamiltonian of section IIIB in [l|) 
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Here 7 is known as the Barbero-Immirzi parameter, a real valued ambiguity parameter of LQG than can be constrained 
from black-hole entropy considerations. The parameters 6b, Sc are functions of the triad variables, the exact form given 
in [l|. These parameters directly measure the magnitude of the loop quantum effects. In the limit of 6b, 6c -^ 0, the 
classical behavior is recovered. 

Choosing the lapse function to be A^ = •j y/p^6b / {sin b6b) , the complete set of equations derived from the effective 
Hamiltonian ([6]) are given by 

sin c6c 1 sin 6bb 

2 6r 



Pb 



6c 

2 6c 
2pc cos c6c 

1 ( sin 
^2 
1 
2 



c cos c6c 



b cos b6b I 1 



J'S! ^ , j'6l 1 



2 6c 



6b 
cos b6b { 1 



1^6b 



sin b6h 



sin c6c 



-^ c cos c6c 

Ob 



Pb- 



(7) 
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With this choice of lapse, the classical singularity would correspond to T — > —00. Furthermore, the system is a 
constrained type whereby the Hamiltonian should vanish, providing an additional equation of motion. 

The results of a numerical solution to the equations of motion are shown in figures [T] The results indicate that 
Pc tends to a constant value as T tends to minus infinity while pb grows as an exponential as in the Nariai case, see 
Eqs. dS]). This is indicative of the space-time metric approaching a Nariai type universe as the classical singularity is 
approached. 
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FIG. 1: In the left figure Pc{T) is plotted which tends towards a constant as T ^ —00. The right figure shows Pb{T) which 
grows like exp(— aT) asymptotically. This asymptotic behavior is indicative of a Nariai type metric. 



In the limit oi T —f —00, there is an exact asymptotic solution to the quantum equations of motion to which the 
solution shown in figure [T] tends to. This solution is characterized by 



Pb^PbS 
Pc = Pc, 



(11) 



where the barred quantities and a are constants. It is clear that this asymptotic solution corresponds to a Nariai type 
solution ([5]) in the large negative T limit. From the equations of motion we can determine the asymptotic values of 
the constants pc, b, a, and c/pb- A root find algorithm can be used to determine the values of the constants which are 
given by [l| 

b w 0.156, pc w 0.182 /li, aw 0.670 

c/pb~ -2.290 mil, TV « 0.689, (12) 

where N is the asymptotic value of the lapse which also behaves as a constant. It is noteworthy that the fixed 
two-sphere radius is Planckian in length since the value of pc approaches a constant Planckian value of 0.182/pj. 
Furthermore this value is independent of the mass of the black hole. 

Given that the classical Nariai solution is unstable to perturbations of the two-sphere radius, we can now ask 
the same question in the context of the loop quantum equations. We consider perturbations around the asymptotic 
solution (fTTj) of the form 

c = ce-"^(l-f£c(i)(r)), (13) 

Pc=Pc + ep['HT), (14) 

b = b + eb^^\T), (15) 

Pb=pbe-^^{l+epi'\T)). (16) 

Now, we insert these expressions into the equations of motion and linearize the equations with respect to e. Since 
we perturb about the Nariai background the zero order equations are identically satisfied. Since the four equations 
of motion are constrained due to the vanishing of the Hamiltonian, we are effectively left with three independent 
equations. This linear system of differential equations describing the perturbations then becomes 





-6.4778 2.2408 \ [c^^^^ 
1.5708 -0.5717 p'^^ 
16.0296 J L(i) 



(17) 



and is solved by 



c(i) = (4.8426 - 3.84266'^'^ cos(wT) -I- 2.5552e'''^ sin(wr))4^^)C, (18) 

p(i) = (cos(wT) - 0.6271 sm{LuT))e^^^, (19) 

pW = (4.8426 - 3.8426e*^ cos(wT) - 4.8510e''^ sin(wr))c[,^^)^, (20) 

where w = 3.0390, S = 0.3350 and where ^ denotes the initial perturbation. Note that we are interested in the 
dynamical behavior of the perturbation close to the classical singularity, located at T = — oo. Therefore, one can 
immediately conclude that the perturbations damp off. A plot of the perturbation behavior is shown in figure [2l 
showing the damped oscillatory behavior for T — > — oo. 

We have thus shown that the Nariai behavior near the classical singularity is a stable solution. The quantum 
interior will not decay into a separate Schwarzschild-de Sitter universe in contrast to the the classical Nariai universe. 
However, the most important issue of how the loop modifications affect a real inhomogeneous collapse scenario remains 
an open problem. 
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FIG. 2: In the left figure we plotted the perturbation 5pc as a function of T. We see that the amplitude of the initial perturbation 
S = 0.01 decreases towards the location of the classical singularity and eventually vanishes. The right figure show the behavior 
of the perturbation Sc which also shows a damped oscillation. 
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